General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



JSC-09379 


NASA TECHNICAL MEMORANDUM 


r. 




\ 



NASA TM X-S8193 
May 1976 


THE ADAMS FORMULAS FOR NUMERICAL INTEGRATION 
OF DIFFERENTIAL EQUATIONS FROM 1ST TO 20TH ORDER 


(HASi-TH-X-58182) THE ADIHS FOBMOLAS FCB 
NUKEBICAL ISTEGBATICB OF EIFEBBENTIAL 
ECO*"IOHS FBCn 1ST TO 20TH CBDEB (MASA) 

37 p HC $4.00 CSCL 12A 


B76-26919 


Unclas 

42323 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
LYNDON B JOHNSON SPACE CENTER 
HOUSTON, TEXAS 77058 



2. Government Accei«ion No. 


3. Reciptent't CaUlog No. 


1. Report No. 

NASA TM X- 58182 

4. T itio end Subtitle 

THE ADAMS FORMULAS FOR NUMERICAL INTEGRATION 
OF DIFFERENTIAL EQUATIONS FROM 1ST TO 20TH ORDER 

7. Author(t) 

James C . Kirkpatrick 

9. Performing Orgmintion Neme end Addreti 

Lyndon B . Johnson Space Center 
Houston, Texas 77058 


12. Sponsoring Agency Name and AddreM 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


15. Supplementary Notes 


5. Report Date 

May 1976 

6. Performing Organization Code 

JSC-09379 

8. Performing Organization Report No. 

10. Work Unit No. 

986-16-00-00-72 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Memorandum 

14. Sponsoring Agency Code 


16. AMract 

The Adams-Bashforth predictor coefficients and the Adams-Moulton corrector coefficients for the 
integration of differential equations are presented for methods of 1st to 20th order . The order of 
the method as presented in this report refers to the highest order difference formula used in 
Newton's backward difference interpolation formula, on which the Adams method is based. The 
Adams method is a polynomial approximation method derived from Newton's backward difference 
interpolation formula . The Newton formula is derived and expanded to 20th order . The Adams 
predictor and corrector formulas are derived and expressed in terms of differences of the deriv- 
atives, as well as in terms of the derivatives themselves. All coefficients are given to 18 signifi- 
cant digfits . For the difference formula only , the ratio coefficients are given to 10th order . 


17. Key WordJ (Suggested by AuthorlsM 
Predictors 

Formulas (mathematics) 

Coefficients 

Polynomials 

Measure and integration 


Runge-Kutta method 
Interpolation 
Extrapolation 
Differential equations 


18. Distribution Statement 

STAR Subject Category: 
64 (Numerical Analysis) 


19. Security Oanif. (of this report) 

20. Security Classif. (of this pugel 

21. No. of Page* 

22 . Price* 

Unclassified 

Unclassified 

37 

$4.00 


JSC Form 1424 (Rev Nov 75) 


‘For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA — JSC 

























NASA TM X-58182 


THE ADAMS FORMULAS FOR NUMERICAL INTEGRATION 
OF DIFFERENTIAL EQUATIONS FROM 1ST TO 20TH ORDER 


James C. Kirkpatrick 
Lyndon B . Johnson Space Center 
Houston, Texas 77058 


CONTENTS 


Section Page 

SUMMARY 1 

INTRODUCTION 1 

ANALYSIS 2 

CONCLUDING REMARKS 9 

REFERENCE 9 

APPENDIX A — NEWTON'S BACKWARD DIFFERENCE FORMULA 

EXPANDED TO 20TH ORDER A-1 

APPENDIX B — ADAMS PREDICTOR AND CORRECTOR 

FORMULAS FOR 1ST- TO 20TH-ORDER SYSTEMS 

IN TERMS OF DERIVATIVE VALUES B-1 

APPENDIX C — ADAMS PREDICTOR AND CORRECTOR DIFFERENCE 

FORMULAS C-1 

APPENDIX D — EXPANSION OF DIFFERENCE FORMULAS IN TERMS 

OF TABULATED VALUES D-1 


i'ui' CilDING PAGE BLANK NOT FILMED 


V 


THE ADAMS FORMULAS FOR NUMERICAL INTEGRATION 


OF DIFFERENTIAL EQUATIONS FROM 1ST TO 20TH ORDER 

By James C. Kirkpatrick 
Lyndon B . Johnson Space Center 


SUMMARY 


The Adams integration method is a fast , accurate polynomial approximation 
method for integrating differential equations numerically . The method , which was 
derived from Newton’s backward difference interpolation formula, requires only two 
derivative evaluations for each integration step . The formulas are derived in terms 
of differences of the derivatives , as well as in terms of the derivatives themselves . 
All coefficients are given to 18 significant digits. For the difference formula only, 
the ratio coefficients are given to 10th order . The binomial and factorial coefficients 
associated with the difference formulas are also given to 20th order. 

The accuracy of the Adams integration method is a function of (1) the accuracy 
of the first n + 1 derivative values, which must be initially supplied, and (2) the 
order of the method chosen in relationship to the word length of the computer . (The 
higher the order of the method used , the smaller the size of the integration step that 
can be used; however, an order that is too high or too low can lead to serious errors 
and cause the solution to diverge . ) 


INTRODUCTION 


The Adams formulas constitute one of the fastest and most accurate methods for 
numerical integration of differential equations; however, there are two major draw- 
backs: the formulas are not self- starting, and the amount of labor associated with 
the computation of the higher order formulas is nearly prohibitive. The first draw- 
back can be overcome by using a self-starting method to initiate the solution. Of the 
self-starting methods, the fourth-order Runge-Kutta method (with the Gill modifica- 
tion) and the Fehlberg 7/8 method are the most common. The second drawback can 
be overcome by judicious programing. It is hoped that the results presented in this 
report will eliminate the need for programing the solution for the computation of re- 
quired coefficients . 

The Adams method is a polynomial approximation method derived from New- 
ton's backward difference interpolation formula. The Newton formula is derived in 
the analysis section, and the results are expanded to 20th order in appendix A. The 
Adams formulas are also derived in the analysis section , and those results are ex- 
panded for use in 1st- to 20th-order systems in appendix B . Appendixes C and D 
provide, respectively, the expansions of the Newton formula used in the derivation 


of the Adams formulas and the difference formulas in terms of the tabulated function 
values . Both expansions are carried out to 20th order . 


ANALYSIS 


Newton's backward difference formula is used for interpolating a value of a 
function y = f(x) near the end of an ordered set of tabular values y. . Thus, if 

y = f(x) denotes a function that takes the values y^^, yj^.j. yj^_2 Yj for the 

n + 1 equidistant values yn_2 the independent variable x, 

then finding a suitable polynomial F(x) to replace the tabulated function over a 
given Interval can be accomplished (ref. 1) by 


n-1 


n 


= «o *2 TT (* ■ *k) 


i=0 k=n-i 


= »0 * (* - *n) * *2 (* ■ *n-l) (* - *n) 

^»3(*-V2)(*- Vl)(*-*n) 

* ®r(* • *n-r+l)(* ' *n-r^-2) ‘ ' *n) 




( 1 ) 


The coefficients a^ , a^ , 82 , . . • * will be determined to make F (x ) = y . 

=yn-r =V2 =Yr 

Substituting successive values x^, x^_^ x^ for x, together with the 

corresponding tabular function values y^, Yj^f in equation (1) gives 


2 


( 3 ) 


■'n-1 0 1 ^ n-1 n) 

■'ll 1^ n-1 nj 

Because the are equidistant , = h and 

. . (yn ■ yn-l) _ 

“l - ^ h R- 


Similarly , 

yp-2 = *0 * “l(*n-2 ■ *n)* ®2(*n-2 ' *n-lX*n-l ' *n) 
■ ^n-l 

= y„ + h ~ ■*■ ®2 

= -^n * =^n-l * 

„ _ yn - %-l ^ yn-2 _ Vn 
* 2ti^ 2h^ 

Equation (5) may be verified by recalling that 

^2yn = '^lyn - '^iyn-1 

= (yn-yn-l)-(Vl-yn-2) 

= yn'*yn-l*yn-2 


> (5) 


( 6 ) 


A list of the higher order differences in terms of the tabulated function values is 
given in appendix D . 
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yn-3 = *0 * »1 ( V3 - *n) * »2 (*n-3 ' *n-l)(*n-3 ' *n) 

* *3(*n-3 - *n-2)(*„-3 ‘ *n-l)(*n-3 ' *n) 

= <-’•>) * Vz j j 

2h 

+ ag(-h)(-2h)(-3h) 

= y^a - 3 t 3) + y^_j(3 - 6) + yn.j(3) - Sh’aj 

= yn-®Vl*»yn-2-®‘>’«3 
a = ^n-1 * ^^11-2 ~ ^11-3 _ ^3^11 


In general , for 1< i < n , 


A.y 
■ ilh* 


Equation (1) can now be written as follows. 


tr(^) 


k=n-i 




Vn /* ■ - *2 


) ■ ■ -(^) 
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Let 


u = 


* ' 


h 


or X - A + hu 
n 


( 10 ) 


Then, because - h, „ = x_ ~ 2h, .... the above factors can be written 

as follows. "'1 " "■* " 


X 


n-1 


K 




+P-U+1 


X - 





u T 2 




> ( 11 ) 


X ■ Xi _ X - 

"T H 


u + n - 1 


/ 


Substituting equation (11) in equation (9) gives 


" A,v >'-1 


F(x)=y„*2 4f TT 

k=l i=0 

^ u(u + 1) (u + 2) 


3! 


"^ 3^0 * • • • 


, u(u + l)(u + 2) . . . (a + n 

* SI — 


- 1) 


A y 
n'^n 


( 12 ) 


Equation (12) is Newton's backward difference formula. It has been expanded to 
20th order in appendix C . The formula can be placed in terms of the function values 
if the difference formulas in appendix D are used . 

If the first n 1 derivative values of the derivatives of a differential equation 
are known at n 1 consecutive , equally spaced values of the independent variable, 
then a polynomial can be fitted to these data and the equation can be integrated by 
integrating the polynomial. (The appropriate polynomial for this purpose is that 
given by equation (12) ) . Integration can be done because, in the numerical integra- 
tion of differential equations , the solution is always begun from a set of given values 
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(initial conditions) and constructed from that point forward . Hence . the values of 
the function immediately behind at a given point of the solution are always known • 
but the values ahead are unknown . The problem of finding the next value ahead is 
solved by the Adams predictor and corrector formulas , which are derived as follows. 

By applying Newton's formula to a set of ordered values of the derivatives of 
a function to be integrated (e.g. , writing y' dy/dx for y in equation (12) and 
using the expanded terms from appendix C) , 


y' = y^^ + u A^y^ + ^u^ + u^ + ^u^ + 3u^ + 2u^ ■ 

+ (u* + 6u^ + llu^ + + . . . 


(13) 


Because the change in y for any interval is given by the formula 


>k+l 


Ay = 


dx 


-I VO. 

*k * 


?k+l 


y'dx 


(14) 


equation (13) can be used to compute the change in y over any interval wherein 
dy/dx is continuous. Therefore, over any interval Xj^^^ - Xj^, 


Xy ^ J l^yj^ + u A^y^ +^u^ + u) 




/ 3 2 

+ ^u + 3u + 2u 




(u^ + 6u^ + Uu^ + 6u^— ^ + . . .j 


dx 


(15) 


Because 


x = x^ + hu 
dx = h du 


( 16 ) 
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equation (15) becomes 


[ 2 / 3 2\A„y' 

VA^It *t)-P 

J 11 


after integration with the change in limits 


*k-H ■* “k*! 


and 


In the inter- 



CIS) 


Integrating equation (17) between the limits 0 and 1 g^ves the following predictor 
formula . 


m ^S^n ^ TW ^7^n 3525500 ^S^n 

. 25713 * . 26842253 * . 1 

^ 55555 Vn ^ 55550525 "^lO^n ^ • • • I 


(19) 


In the interval x , - x , 
n-i n 


X - X 

_ n n 

0 

'^k+l R 

■ R " 

_ n-1 n 

-h 

'he h 

" TT 


= -1 


( 20 ) 
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Integrating equation (17) between the limits *1 and 0 gives the following corrector 
formula. 


= •> [vn - J - 0 - n AjyA - m 

■ iro ■ ram ' nrai ^?yn ‘ rairara ^jyn 

«183 . 32504331 , 1 

■ irarara ^o^n ' mraisra '^loyn ' ■ • J ^ ' 

Equations ( 19 ) and ( 21 ) are the predictor and corrector formulas of the Adams 
method and are to be used as follows . Use equation ( 19 ) to compute a predicted 
value , 


y 


n+1 



( 22 ) 


Place a bar over to indicate that it is a first approximation obtained by extrap- 
olation. Use this value of y^^j to compute the value of the derivative With 

the value of the y^^j^ derivative thus obtained, compute a corrected value of 

by using equation (21) with the n and n-1 replaced by k and k-1, wherein 
k = n + 1. 


= Yn + h \ ^l^ii+l ■ U ^2^0+1 ‘ W ^3^n+l ’ ' ' ] 

Equations (22) and (23) have been placed in an easy-to-use form in appendix B. 
The equations given in appendix B are written in terms of the derivative values at 
the various intervals ranging from a 1st- to 20th-order system. For an nth-order 
system .n + 1 derivative values are required . 
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CONCLUDING REMARKS 


The Adams integration method is a fast, accurate method for integrating dif- 
ferential equations numerically. It requires only two derivative evaluations for each 
Integration ste]^, hut addition^ time and storage must be dedicated to continually up' 
dating or reordering the table of past derivative values. When large sets of differen- 
tial equations are integrated by using a high-order mediod . the amount of time re- > 

' quired to update the derivative table can become appreciable. When this occurs, 
the Adams formula is preferred for its speed. 

The accuracy of the Adams integration method is a function of (1) the accuracy 
of the first n + 1 derivative values, which mustbe initially supplied, and (2) the 
order of the method chosen in relation to the word length of the computer. The high- 
er the order of the method used, the smaller the size of the integration step that can 
be used; however , an order that is too high or too low can lead to serious errors 
and cause the solution to diverge . For single-precision work , no order higher than 
seventh or eighth should be used with a nine-digit word length . For double- 
precision work, 13th is probably the highest order that should be attempted with a 
20-digit word length. A lOth-order method can be used successfully witti a 12-dl^t 
word length, and an order up to 20th can be used successfully with a 24-digit word 
' length. 

Because the data in this report were computed using floating-point arithmetic, 
only the first 16 digits of each 18-digit number are recommended without reservation. 


Lyndon B. Johnson Space Center 

National Aeronautics and Space Administration 
Houston, Texas, April 20, 1976 
986-16-00-00-72 
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APPENDIX A 


NEWTON'S BACKWARD DIFFERENCE FORMULA EXPANDED TO 20TH ORDER 
The expansion of Newton's backward difference formula to 20th order is as 


follows . 



k-1 

. . 4. V 

y -yA*2-( -et 

(u + i) = y^^ + uAjy^^ 

k=l 

=0 

/ 2 

+ iu +u; 2 



+ (u’ + 3u* * 2u) ^ 

* (u< * 6u’ + nu^ + 6u) ^ 


+ (u® + lOu^ + 35u® + 50u^ + 24u) 

A y' 

+ (u® 15u® + 85u^ + 225u^ + 274u^ + 120u) 

+ (u*^ + 21u® + 175u® + 735u^ + 1624u^ + 1764u^ + 720u) 

A y * 

+ (u® + 28u^ + 322u® + 1960U® + 6769u^ + 13132u® + 13068u^ + 5040u) 

+ (u® + 36u® + 546u^ + 4536u® + 22449u® + 67284u^ + 118124u® + 109584u^ 

" “osaou) 

+ (u^® + 45u® + 870u® + 9450u^ + 63273u® + 269325u® + 723680u^ + 1172700u‘ 
+ 1026576u^ + 362680u)g^^j^ 


A-1 


+ (u^^ + 55u^® + 1320U® + 18150U® + 157773u^ + 902055u® + 3416930u® + 8409500u^ 


+ 127535760u^ 10628640u^ * 3628800u) 

+ + 66u^^ + 1925u^° + 32670U® + 357423u® + 2637558u^ + 13339535u® 

A y* 

45995730u^ + 105258076u^ + 150917965u^ 120543840u^ + 39916800u) 4 7 ^ 001 ^ 

+ (u^® + 78u^^ + 2717u^^ + 55770u^° + 749463u® + 6926634u® + 44990231u^ 

-I- 206O701S0u^ -t- 657206836u^ + 1414014888u^ + 1931559552u^ + 1486442880u^ 

+ 477J')1600u) g2a70206flg 

+ + 91u^^ + 3731u^^ + 91091u^^ + 1474473u^° + 16669653u® + 135036473u® 

+ 790943153u^ + 3336118786u® + 9957703756u® + 20313753096u^ + 26596717056u® 

2 \ ^ 14 ^^ 

+ 19802759040u^ + 6227020800u) 

+ (u^® + 105u^^ + 5005u^® + 143325u^^ + 2749747u^^ + 37312275u^® + 368411615u® 

+ 2681453775U® - 1 - 14409322928u^ + 56663366760u® + 159721605680u® + 310989260400u^ 

3 2 \ ^IS^n 

+ 392156797824u'^ -i- 283465647360u^ -t- 87178291200u ) 15076 ^ 43^^0 

+ (u^® + 120u^® + 6580u^^ + 218400u^® + 4899622u^^ + 78558480u^^ + 928095740u^® 

+ 8207628000U® + 54631129553u® + 272803210680u^ + 1009672107080u® 

+ 2706813345600U® + 5056995703824u^ + 6165817614720u® + 4339163001600u^ 

+ 1307674368000u) gOffgg jg ftggfiOO - O 


A-2 


+ + 136u^® + 8500u^® + 323680u^^ + 8394022u^® + 156952432u^^ + 2185031420u^^ 

+ 23057159840u^° + 185953177553u® + 1146901283528u® + 5374523477960u^ 

+ 18861567058880u^ 48366009233424u^ + 87077748875904u^ + 102992244837120u^ 

+ 70734282393600u^ + 20922789888000u) SSSest ' fifflTOOO 

+ + 10812u^® + 468180u^® + 13896582u^^ + 299650806u^® + 4853222764u^^ 

+ 60202693980u^^ + 577924894833u^® + 4308105301929u® + 24871845297936u® 

+ 110228466184200u^ + 369012649234384u^ + 909299905844112u^ + 1583313975727488u^ 

A y* 

+ 1821602444624640u^ + 1223405590579200u^ + 355687428096000u) 6402373705728606 

+ (u^® + 171u^® + 13566u^^ + 662796u^® + 22323822u^® + 549789282u^^ + 10246937272u^^ 

19 11 in 9 

+ 147560703732U + 1661573386473u + 14710753408923u + 102417740732658u 

•t- 557921681547048u^ + 23S3125040549984u'^ + 7551527592063024u^ 

+ 17950712280921504U® + 30321254007719424u^ + 34012249593822720u^ 

+ 22376988058521600u^ + 6402373705728000u) 12164glt64666 ' 3306 ' 6 

+ (u^° + 190u^® + 16815u^® + 920550u^'^ + 34916946u^® + 973941900u^^ + 20692933630u^^ 
+ 342252511900u^^ + 4465226757381u^^ + 46280647751910u^^ + 381922055502195u^° 

+ 2503858755467550u^ + 12953636989943896u^ + 52260903362512720u'^ 

+ 161429736530118960U® + 371384787345228000u® + 610116075740491775u^ 

-t- 668609730341153279u^ + 431565146817638400u^ 

+ 121645100408832000u) 2432902008176640000 
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APPENDIX B 


ADAMS PREDICTOR AND CORRECTOR FORMULAS TOR 
1ST- TO 20TH-ORDER SYSTEMS IN TERMS OF DERIVATIVE VALUES 

The predictor and corrector methods using derivative values directly are as 
follows . 


IST-ORDER SYSTEM 
Predictor 

Corrector 

2ND-ORDER SYSTEM 
Predictor 

^n+1 ^ ^n AT(0.416666666666666667y^_2 - 0.133333333333333333 + 001y^_^ 
+ 0.191666666666666667 -* OOly^^) 

Corrector 

^n+1 "" ^n AT(-0. 833333333333333335 - 001yj^_^ + 0.666666666666666667y^^ 

+ 0.416666666666666668yj^^^) 


B-1 


3RD-0RDER SYSTEM 


Predictor 

^n+1 ^T(-0.375yj^_3 + 0.154166666666666666 + OOly^^.j 

- 0.245833333333333333 OOly' , •«- 0.229166666666666666 + OOly’ ) 

n“i n/ 

Corrector 

^n+l ^ AT(0. 416666666666665 - COly^^.j ‘ 0.208333333333333333yj^_^ 
+ 0.791666666666666666y^ + 0.375y^^^) 


4TH-ORDER SYSTEM 


Predictor 

^n+1 " AT(0.348611illlllllllllyj^_4 - 0.176944444444444444 + 001y^^_3 
+ 0.363333333333333333 + OOly’ „ - 0.385277777777777777 + OOly' , 

n~4 n“i 

+ 0.264027777777777777 + OOly^^) 

Corrector 

^n+1 "" AT(-0. 263888888888888892 - 001yj^_g + 0.147222222222222224y^_2 

- 0.366666666666666669y^_^ + 0.897222222222222223y^ 

+ 0.348611111111111113yj^^^) 
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5TH-ORDER SYSTEM 


Predictor 

^n+l "" AT(-0.329861111111111111yj^_5 + 0.199791666666666666 + 001yj^_4 

- 0.506805555555555554 -f 001y'_„ 0.693194444444444443 •<- OOly* » 

n u n A 

- 0.550208333333333331 + 001yj^_^ + 0.297013888888888888 + OOly^^) 

Corrector 

yn^l =yn + AT(0.1875 - 001y^_^ - 0.120138888888888889yj^_g 

+ 0.334722222222222223y^^_2 - 0.554166666666666668y^_^ 

+ 0.990972222222222222y^ + 0.329861111111111114yj^^j) 

6TH-ORDER SYSTEM 
Predictor 

^n+1 ^ AT(0.315591931216931216yj^_g - 0.222341269841269840 + 001yj^_5 

+ 0.673179563492063490 + 001y^_^ - 0.113798941798941799 + 002yj^_3 

+ 0.116658234126984127 + 002y' „ ~ 0.739563492063492062 + OOly' , 

n~2 n-i 

+ 0.328573082010582010 + OOly^) 

Corrector 

^n+1 ^ AT(-0. 142691798941798941 - 001y^_g + 0. 104365079365079364y^_^ 

- 0.334176587301587300yj^_3 + 0.620105820105820105y^_2 

- 0.768204365079365079y^_^ + 0.107658730158730159 + OOly^ 

+ 0.315591931216931219y^^^) 
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7TH-0RDER SYSTEM 


Predictor 

^n+l "" AT(-0.304224537037037036y^_,y + 0.244516369047619047 + OOly^^.g 

- 0.861212797619047617 + 001y^_g + 0.173796544312169311 + 002y^_^ 

- 0.220277529761904761 + 002y' „ + 0.180545386904761904 002y' „ 

''n~3 ■'n-^ 

- 0.952520667989417986 001y^_^ + 0.358995535714285713 + OOly^^^ 

Corrector 

^n+1 ^ AT(0. 113673941798941799 - 001y^_g - 0.938409391534391534 - OOly^.g 

+ 0.343080357142857143y^_4 " 0.732035383597883596y^_g 

-<-0.101796461640211640 -<- OOly* » - 0.100691964285714286 -<- OOly' , 

n-<6 •'n-i 

-I- 0.115615906084656084 -<■ OOly^^ -<- 0.304224537037037040yj^^^) 

8TH-ORDER SYSTEM 
Predictor 

^n-<-l ^ AT(0.294868000440917107y^_g - 0.266316854056437389 -<- 001y^_^ 

-<■0.107014677028218694 -<- 002y^_g - 0.251247360008818341 -<- 002y^_g 

-t- 0.380204144620811286 ■<■ 002y^_^ - 0.385403610008818341 -<^ 002yj^_3 

-<-0.263108427028218694 -<- 002y' „ - 0.118841506834215167 -<- 002y' , 

n~z n-i 

■<■0.388482335758377424 -<^ OOly^) 
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Corrector 


^n+1 AT(-0. 5653659611992939 - 002y^_^ + 0.862196869488536150 - 001y^_g 

- 0.355823963844797176y^_g + 0.867046406525573189y^_^ 

- 0.138699294532627865 + 001y^_3 + 0.154193066578483245 + 001y^_2 

- 0.126890266754850088 + 001y^_j + 0.123101135361552028 + OOly^^ 

+ 0. 2948680004409 17110y^^j) 

9TH-ORDER SYSTEM 
Predictor 

^n+1 " AT(-0.286975446428571428yj^_g + 0.287764701829805996 + 001y^_g 

- 0.129942846119929453 + 002y^_^ + 0.348074052028218694 + 002y^_g 

- 0.612836422508818339 + 002y^_g + 0.741793207120811284 + 002yj^_^ 

- 0.626462985008818339 + 002y^_g + 0.366419587742504408 + 002yj^_2 

- 0.144669297012786595 + 002y' , + 0.417179880401234566 + OOly' ) 

n-l n/ 

Corrector 

^n+l ^ AT(0. 789255401234567893 - 002yj^ g - 0.803895227072310397 - OOly^!^ ,^ 

+ 0.370351631393298056y’ „ - 0.101879850088183421 + OOly' - 

n^o n“0 

+ 0.186150821208112873 + 001y^_^ - 0.238145475088183419 + 001yj!^_3 

+ 0.220490520282186947 + OOly' 0 . 155303461199294532 + OOly' . 

•'ll- 2 'n-l 

+ 0.130204433972663139 + OOly^ + 0.286975446428571432yj^^^) 


B-5 


lOTH-ORDER SYSTEM 


Predictor 

^n+1 * * ^T(0.280189596443936721yjJ_jQ - 0.308887141086793863 + OOly^^.g 

0.154861788582752124 ••• 002y^_g - 0.466170361852653517 -i- 002y^.^ 

-»■ 0.936472204560485808 -t- 002y^_g - 0.131891420554753887 + 003y^_g 
+ 0.133019135965307840 + 003y^_^ - 0.962690500741542404 + 002yj^_g 
+ 0.492504906142275932 + 002yj^_2 - 0.172688256657180267 + 002y^^_^ 

+ 0.445198840045628238 + OOly^^ 

Corrector 

^n+1 “ AT(-0. 678584998463470693 - 002y^_g + 0.757510538586927482 - OOly^^.g 

- 0.385752772015792851y^_^ + 0.118465362954946289 + OOly^^.g 

- 0.244382699765512266 + 001y^_g + 0.357154240820907487 + 001y^_^ 

- 0.380648324765512264 OOly* . 0.301920720097803430 ^ OOly* » 

n~o n-A 

> 0.185839786130150713 OOly' , + 0.136990283957297846 + OOly' 

n~i n 

+ 0.280189596443936725y^^^) 

IITH-ORDER SYSTEM 
Predictor 

5^n+l ^ AT(-0.274265540031599059yj^_^j + 0.329711053679152637 + OOly^^ ^g 

- 0.181734761126058868 + 002y' . + 0.607399929634890570 •»- 002y' q 

■^n-8 n-o 

- 0.137124664395693041 + 003y^_^ + 0.220357899950647346 + 003y^_g 
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- 0.258602100049352652 •»- 003y*_. + 0.223526764175735529 -»■ 003y'_ . 

n D li’*^ 

<• 0.141522864179368085 + 003y' » + 0.643350953159655414 002y'_, 

H“u n A 

- 0.202857466060656163 + 002yj^_^ + 0.472625394048788144 + OOly^^) 

Corrector 

^n+1 " AT(0. 592405641233766230 - 002y^_jQ - 0.719504705203489922 - 001yj^_g 

+ 0.401574156537264175y’ ^ - 0.136322208005150713 + OOly’ - 

ri“’o I 

+ 0.313959224562089145 + OOly’ - - 0.518074106015512264 + OOly’ - 

ri“0 ii“3 

+ 0.630845647070907485 + 001y^_^ - 0.576142186372655120 + 001y^_3 

+ 0.399667650901374858 001y’_„ - 0.218422096398007856 + OOly' , 

n « n X 

+ 0.143506746010869274 + OOly^^ 0 .274265540031599064y^^j^) 

12TH-0RDER SYSTEM 
Predictor 

^n +1 " AT(0.269028846773648774y^_j2 " 0.350261170131538434 + 001 y^_^j^ 

+ 0.210530144238523454 + 002y^_^Q - 0.773598224028086170 + 002y^_g 

+ 0.193909272116445200 + 003y' 0.350195511040422869 + 003y! „ 

ri“o n I 

+ 0.468940554369498813 + 003y’ „ - 0.471672946694082481 + 003y' - 

n-o n-o 

+ 0.356696043328691671 + 003y' . - 0.200709210469570815 + 003y' , 

n~4 n"o 

+ 0.820909992030263604 + 002y' « - 0.235140927673494016 + 002y' , 

n** « HI 

+ 0.499528278726153021 + OOly^^) 
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corrector 


^n+1 ^ AT(-C. 523669325795028499 - 002yj^.^j + 0.687643755077410822 - OOly^^.^Q 

- 0.417572225545067802y^_g + 0.155364667328632887 + 001y^_g 

- 0.395538524273689820 + OOly^^,^ + 0.728705330591751715 + 001y^_g 

- 0.100194456305011860 -i- 002y^_g -i- 0.104559175310057006 002y^_^ 

- 0.835358502641194226 + OOly^.g -t- 0.514874902576281128 001yj^_2 

- 0.252984271900479736 OOly' , + 0.149790777020409616 OOly' 

n~i n 

+ 0.269028846773648779y^^j) 

13TH-0RDER SYSTEM 
Predictor 

^n+1 ‘ AT(-0.264351348366606509yj^_j3 + 0.370559637553953339 + OOly^^.^j 

- 0.241220168739106920 + + 0.968575000587018070 + 002y^^_jQ 

- 0.266371036484932270 003y* . 4^ 0.534129457464267776 003y' . 

n-w n~8 

- 0.803822424837519638 + 003y^.^ 0.922567468166595582 4- 003y^,g 

- 0.811883132041905058 4> 003y^.g 4- 0.545707257410815325 4> 003y^.^ 

- 0.276313696102420276 4- OOOy^^.g + 0.102710404375621668 4- 003yj^_2 

- 0.269506602961152862 4- 002y^_^ 4^ 0.525963413562813672 4^ 00 ly^) 

Corrector 

^n4-l ^ AT( 0.487749840704226446 - OGOy^^.^g - 0.660441725494997230 - 001y^_j^ 

4- 0.433609251257037710y^_^g - 0.175533676995915544 4- 001y^_g 
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0.489805803432154596 ^ OOly^^g - 0.997532539260029256 OOly^.^ 

+ 0.153136405724020429 * 0O2y^_g - 0.180460328969357118 + 002y^_g 
^ 0.164758579808690949 002y^_^ - 0.116979963874471613 + 002y^,^ 

+ 0.648651357017689891 + 001yj^_2 ' 0.289468759475409399 + 001y^_^ 

+ 0.155871525849564560 + OOly^^ + 0.264351348366606514y^^^) 

14TH-0RDER SYSTEM 
Predictor 

^n+1 "^n ^T(0.2601.'»6396127601036y^_j4 - 0.390626089415302101 + OOly^^^^g 

+ 0.273780084231512276 + 002yj^ ^2 " 0.118811565064357469 + OOSy^^,^^ 

+ 0.357054032580430444 + 003yj^_^p - 0.787164101532389544 + 003yj^_g 

+ 0.131531905503545369 + 004y' „ - 0.169661053634744639 + 004y’_, 

n”o n I 

+ 0.170375706573778149 + 004y' - - 0.133268619708936233 + 004y' , 

n~D n “0 

+ 0.806103789934543962 + 003y^_^ - 0.371003344292867054 + 

+ 0.126382816423233362 + 003y’ „ - 0.305925698419017007 + 002y’_ , 

+ 0.551977053175573776 + OQly^) 

Corrector 

^n+1 AT(-0.421495223900547286 - 002y^_jg + 0.636868297531188845 - 001yj^_j2 

- 0.449604826298997753yj^_jj + 0.196785186625502983 + OOly^-io 

- 0.597450396120363375 + 001yj^_g + 0.133363924168105026 + 002yj^_g 

- 0.226328272663337275 + 002y' „ + 0.297793566566688258 + 002y' - 

n-< n-D 
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- 0.307035344707191487 + OOZy^^g + 0.249141023033580516 + 002y^.^ 

- 0.159171635786916397 + 002yjJ_3 + 0.802075618517489103 + OOly^J^g 

- 0.327824824850359202 + 001y^_^ + 0.161772458984172222 + OOly^ 

+ 0.260136396127601041yjJ^j) 

15TH-ORDER SYSTEM 
Predictor 

^n+1 " ^ AT("®-256309496574389152yjj_jg + 0.410477884474343831 + OOly^^^^ 

- 0.308187580344638819 + 002y^_^3 -f 0.143998829364498292 + 003y^_^2 

- 0.468674127888398661 'fOOSy* + 0.112675145079332106 + 004y' 

•'n-ii 'n-10 

- 0.206989313188720724 * 004yl ^ + 0.296467066549164788 + 004y' 

n“» n-p 

- 0 . 334596214680364058 + 004y^_^ + 0 . 298658609609259920 + 004y^_^ 

- 0.210238361530225295 -i* 004y^^g -f 0.115596625275858515 -i- 004y^_^ 

- 0.487624165234214118 + OOSyJ^ ^ + 0.153295313563544223 -t- 003y* « 

n—o ii“£ 

~ 0.344372122905175380 + 002y^_^ + 0.577608002833012691 + OOly^^ 

Corrector 

^n+l " AT(0.382689955321188443 - OOZy^^^j^ - 0.616184455371837393 - OOly^J.^j 

+ 0.465511282840366749yj^_j2 - 0.219084412301040517 + OOly^.jj 

fk ■ 

-i- 0.719156975638925207 -I- OOly^.^Q - 0.174666833194989227 + 002y^.g 

+ 0.324900246806359841 -i- 002y^_g - 0.472589258912522038 002y^^^ 

-I- 0.544054552815873020 + 002y^_g > 0.498571667345446283 > 002y^_g 
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+ 0 .364063717216533405 + 002y* . * 0 . 211408814688258610 + 002y* » 

‘n-a 

+ 0 .976109548188629844;** ^ 0.368007270159083988 + 001y^_^ 

+ 0.167512808313990048 ^ OOly^ 0.256309496574389157y^^^) 

16TH-ORDER SYSTEM 
Predictor 

y^^j = y^^ + AT^O . 252812146729039235yj^_^g - 0 . 430130384423901690 + 001y^_ jg 
+ 0.344422364522281464 + - 0.172393560202725853 + 003yj^_jg 

+ 0.601116936411349698 + 003yj^_j2 " 0.157295758480084203 + 004yj^^j^j^ 

+ 0.315127112179946725 + 004y^_j^Q - 0.496216409046741608 + 004y^_g 
+ 0.621836299389438282 *f 004y^,g - 0.623813310538384942 + 004y^.^ 

•► 0.501110576709874538 + 004y^,g - 0.3206667Q7221469633 •► 004y^_g 
•► 0. 161608435980543656 •► 004jr^_^ - 0.629198967402476089 003y^_3 

+ 0.183632771171028931 + 003yj^_2 " 0.384822066381821657 + 002yj^_j 

i 

•► 0.602889217505916614 + OOlyj^) 

Corrector 

^n+1 " * ^^(-0.349734984534991753 - 002y^.^5 + 0 .597844970788105648 - OOly^^^j^ 

- 0.481300426979173843y^.j3 + 0.242402719623632056 + 001y^_^2 

- 0.855602084154725506 001y^_jj + 0.224679938808776918 002y^_jg 

- 0.454734608810610621 + 002yjj^_3 0.724997069114390405 + 002y^_g 

- 0.922698184009056423 + 002yjJ^^ + 0.944151375123903584 •► 002yj^_g 
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- 0.778639442961067678 + 002y^^_g + 0.516827958461417802 + 002yj^_^ 

- 0.275060581873627117 + 002y^_g + 0.117205113952822522 + 002yj^_2 

- 0.409975468303282999 + OOly' , + 0.173108568066549916 -h OOly’ 

n“i n 

+ 0.252812146729039240y^^^) 

17TH-0RDER SYSTEM 
Predictor 

^n+1 AT (-0. 2495976502977 15667yj^_^^ + 0.449597220179020557 + 001y^_^g 

- 0.382465842847283475 + 002yj^_^g + 0.204168638654674800 + OOOy^^.^^ 

- 0.766435967911289140 + 003yj^_^3 + 0.214862719645361424 + 004y^^_^2 

- 0.466197810488537112 + 004yj^_^j + 0.800544622478944153 + 004y^^_^g 

- 0.110298829692048839 + 005y’ . + 0.122860818726318507 ^ 005y* . 

n~y n“o 

- 0.110923082083738237 + 005y' _ + 0.810012628718327447 + 004y' « 

n-7 *^n-b 

- 0.475117733225696087 + 004y' - + 0.221012676751399984 + 004y' . 

n~& n~4 

- 0.798925369604922742 + 003y' - + 0.217578051611518262 + 003y' „ 

n-3 'n-Z 

- 0.427253666932433320 + 002y' , + 0.627848982535688180 + OOly* \ 

n-i n/ 

Corrector 

y^+1 = y^^ + AT(0.321449643132356720 - 002y^_^g - 0.5814378917785055?,^ - 001yj^_^g 
+ 0.496956011738815704y^_j^ - 0.266715800027919954 + 001yj^_^3 
+ 0.100745287027864105 + 002y^_^2 ' 0.284473247585774889 + 002yj^_j^ 

+ 0.622506017149331594 + 002y^_^g - 0.107988987477441797 + 003y^_g 
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+ 0.150644115156914959 + 003y' » - 0.170414226646381561 + 003y’ 

n-o n-7 

+ 0.156930664108771093 + 003y' - - 0.117646552130167235 + 003y' - 

n-D •'n-5 

+ 0.715740997631720140 + 002y' . - 0.351565596939128017 + 002y* 

n-4 ’^n-3 

+ 0.139063689685822779 * 002y' » ~ 0.453692619769283512 + OOly' , 

n-2 •'n-1 

+ 0.178573211999799980 + OOly^ + 0.249597650297715674yj^^^) 

18TH-ORDER SYSTEM 
Predictor 

yn+i =yn + AT(0.246628202582257457yj^_18 - 0.468890529677834989 + 001yj^_^^ 

+ 0.422300871968755965 + 002y* 0.239495197591850432 + 003y* 

n-lb *^n-15 

+ 0.958850938556382618 + 003y* 0.287954640763607103 + 004y’ 

■^n-14 ■'n-13 

+ 0.672703314919064167 + 004y’ “ 0.125106740238631324 + 005y* ,, 

n-12 -^n-ll 

+ 0.187974031133838633 + 005y' 0.230209461787542415 + 005y* » 

n~iu •'n-9 

+ 0.230780387612262725 + 005y’ . - 0.189410041273515850 + 005y’ „ 

n-o •'n-7 

+ 0.126785322399203019 + 005y* « - 0.686428777198174276 + 004y' , 

n-D •'n-5 

+ 0.296480906741570766 + 004y' . - 0.100017398291204483 004y' „ 

n-4 *^n-3 

+ 0.255312166606603652 + 003y' „ - 0.471646743397239662 + 002y’ , 

n-J •'n-1 

+ 0.652511802793913926 + OOly^^) 

Corrector 

^n+1 AT (-0.296944771545820966 - 002y^_^^ + 0.566645553095713411 - 001y^_^g 

- 0.512469289642956638yj^_^g + 0.292002534755271478 -t- 
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- 0.117536680095813211 + + 0.355167567288323508 + 002y^^_^2 

- 0.835721521483436930 + 002yj^_^^ + 0.156750305811680223 + 

- 0.237926080610462135 003y^_g + 0.295018663082493113 -»■ 003y^.g 

- 0.300351319779401899 + 003y' „ + 0.251430368205513157 + 003yJ' - 

n~ I n~o 

- 0.172771379519933439 + 003yj^.g + 0.970163277892179544 + 002y^^_4 

- 0.442430697032149232 + 002y^_j + 0.163294383043961770 + 002y^^_2 

- 0.499125169815794119 + OOly' - + 0.183918217887624758 + OOly’ 

n~i n 

+ 0.246628202582257464yj^^j) 

19TH-0RDER SYSTEM 
Predictor 

5^n+l ^T^-0-243872812282820742yj^_jg + 0.488021163595585154 + 001yj^_^g 

- 0.463911561971406966 + 002yj^_^^ + 0.278542842298928894 + OOSy^^.^g 

- 0.118474621800006362 + 004yj^_^g + 0.379460399978102219 + 004y^_^^ 

- 0.949630355049356339 + + 0.190152964144974132 + 

- 0.309430689218232897 + 005y' „ + 0.413258857664462777 + OOby' * 

II" i.u 

- 0.455494288318166560 + 005yj^.g + 0.415104336591864297 + 005y^.g 

- 0.312292673926583565 + 005y^_^ + 0.192952893827777942 + 005y^_g 

- 0.970004083320638234 004y’ . + 0.391006008782392086 + 004y! . 

‘n-5 n~4 

- 0.123648673801409812 + 004yj^_g + 0.297014417506965999 + 003yj^_2 

- 0.517982577730975603 + 002y^_^ + 0.676899084022196000 + OOly^) 
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Corrector 


^n+1 “ AT(0. 275539029943671580 - 002y^^_jg - 0.553218634047558098 - 

+ 0.527836296513249742y^_jg - 0.318244248979713424 + OOly^^.^g 

+ 0.135999181481694252 + 002y^_^^ - 0.437933464114314523 + 002y^_^g 

•H 0.110276006333149324 003y' • 0.222410758556360928 + 003y' 

n-iz n~ii 

+ 0.365008215423706076 + 003yj^_^Q - 0.492463525691827066 + 003y^_g 

+ 0.549556108163858044 + 003y' ^ - 0.508609229391427752 + 003y' „ 

n-o n~7 

+ 0.390268974613530393 + 003yj^_g - 0.247530629124250412 + 003yj!^_5 

•i- 0.129056006191068085 + 003y' . - 0.549229625038316336 + 002y' » 

n-4 ''n-6 

+ 0.189994115045503546 -i- 002y* » - 0.546242343936161959 + OOly* , 

n-A n-i 

+ 0.189153459456554518 + OOly^ + 0.243872812282820749y^^^) 

20TH-ORDER SYSTEM 
Predictor 

^n+1 AT(0.241305789737813504yj^_2Q- 0.506998860703909081 + 001y^_^g 

+ 0.507283116861404171 + 002y^_^g - 0.321479756498248090 + OOOy^^,^^ 

+ 0.144766939357863532 + 004yj^_^g - 0.492595118209512418 + 004yj!^_^g 

+ 0.131476164100186736 + 005y’ , . - 0.282023283709688662 + 005y' 

n-14 n-io 

+ 0.494125867477697802 + 005y^_^2 " 0.714727893661864457 + 005y^_^^ 

+ 0.859085782552457494 + 005yj^ - 0.860791492761798119 + 005y^_g 

+ 0.719077239924587968 + 005y* ^ - 0.499352922131336593 + 005y' - 

n-o n-7 

+ 0.286483017930154456 + 005y' - - 0.134412457973014429 + 005y' , 

•'n-6 "^n-S 
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+ 0.507918663910362727 + 004y^_^ - 0.151157533831520551 + 004y^_3 
+ 0.342862517557150564 + OOSy^^.j - 0.566243735678538303 + 002yj^_^ 
••• 0.701029662995977350 OOly,^) 

Corrector 

^n+1 ^ AT(-0. 256702254500723787 - 002y^_^g + 0.540958411995814733 - 00 

- 0.543056146956131006y^.j7 + 0.345424199782150092 + OOly^^.^g 

- 0.156196667203572017 + 002y^_jg + 0.533990356859616411 + 002yj^.^4 

- 0.143291140255911992 + 003yj^_j3 + 0.309271594022110403 + 003yj^_^2 

- 0.545778588550922682 + 003yj^_^j + 0.796165322083121750 + 

- 0.966736343017184307 + 003yj^_g + 0.980713214823273717 + O03yj!^_g 

- 0.831977059385989507 + + 0.589264562302491472 + O03yj^_g 

- 0.347028422968730952 + 003y' . + 0.168855123728860301 + 003y’ . 

n~D n~e 

- 0.673601867343917011 + 002y^_g + 0.219258172058586058 + 002y^_2 

- 0.595015772291299479 + OOly* , + 0.194287504546568993 -i- OOly* 

n~l n 


n-18 


+ 0.241305789737813511yj^^^) 


APPENDIX C 


ADAMS PREDICTOR AND CORRECTOR DIFFERENCE FORMULAS 


The Adams predictor and corrector difference formulas are expressed as 
follows . 

PREDICTOR DIFFERENCE FORMULA 

y ^1 =y + +0*5 A,y' + 0.416666666666666667 A„yj; + 0.375 A«y' 

•^n+l "^n V n I n 2’'n S'^n 

+ 0.348611111111111111 A.y’ + 0.329861111111111111 A.y’ 

4 n on 

+ 0.315591931216931216 A-y' + 0.304224537037037036 A„y' 

on ( n 

+ 0.294868000440917107 A„y' + 0.286975446428571428 Any' 

on V n 

+ 0.280189596443936721 A,.y’ + 0.274265540031599059 A-,y' 

10 n 11 n 

+ 0.269028846773648774 A^gyn 0.264351348366606509 A^gy^ 

+ 0.260136396127601036 A, .y’ + 0.256309496574389152 A-,y' 

14‘'n 10 n 

+ 0.252812146729039235 A^gy^ + 0.249597650297715667 A^^y^l^ 

+ 0.246628202582257457 + 0.243872812282820742 A^gy^ 

+ 0.241305789737813504 ^20^) 

CORRECTOR DIFFERENCE FORMULA 

^n+1 ^ ^n ^^W+l ' ^l^ii+1 ‘ 0.833333333333333335 - 001 Agyj!^^^ 

- 0.416666666666666665 - 001 A^y^^^ - 0.263888888888888892 - 001 A^y^^^ 

- 0.1875 - 001 Agy^^^j - 0.142691798941798941 - 001 Agy^!^^^ 

- 0.113673941798941799 - 001 A^y^^^^ - 0.935653659611992939 - 002 ^gY^+i 


C-1 


0.789255401234567893 - 002 - 0.878584998463470693 - 002 

0.592405641233766230 - 002 - 0.523669325795028409 - 002 A^2yn^l 

0.487749840704226446 - 002 A^gy^^j - 0.421495223900547286 - 002 A^^y^^^ 
0.382689955321188443 - 002 A^gy^^^ - 0.349734984534991753 - 002 Ajgy^^^^ 
0.321449643132356720 - 002 A^^y^^^^ - 0.296944771545820966 - 002 A^gy^^^ 
0.275539029943671580 - 002 A^gy^^^ - 0.256702254500723787 - 002 Aggy^;^^) 


APPENDIX D 


EXPANSION OF DIFFERENCE FORMULAS IN TERMS OF TABULATED VALUES 


The expansion of the backward difference formula is as follows . 


= K - yA-i 


Vn=yA- 2 yi-i^yn -2 


= K ■ %-i * %-2 ■ yn -3 


= K - %-i * %-2 - %-z * yA-4 


VA = K - %-l * ^K-2 - ^%~Z " %-i - 5A-5 

= K - %-i * ^%-2 - ^K-Z * ^%-i - %-z * K-Z 


= yA ■ ’^A-i * ”yA -2 • ^%-z * ^%-i - ^%-z * %-z - yA-7 


^s^A = yA ■ «yA-i * =*«yA-2 ■ ^%-z * ''K-i - ^%-z * ^K-z - K-i * K-a 


‘^ayk = yk ■ %-i * ^%-2 - ^%-z * ^^%-i - i26yA-5 " ^%-z ‘ ^%-i * %-a 

-yk-a 

^1% = yA - i»yA-i * «yA-2 - ^^%-z * ^^%-i - ^^K-z * ^^%-z - ^^%-i 

* *%-a - ^%-a * yA-io 


D-1 


= yk - * *%-2 - l«5y'_3 * 330y^.^ - 462y'_5 . 462y'_5 - 330y^., 

* l«5y^.j - 55y^_, + lly^.jg - y^_^j 

^125 ^ = yk - * *«Vn-2 ■ ^^K-i * 

* 495y'.g - 220y;.j * 66y^_ig - 12y^.jj . y-.j, 

= yk - ^K-I * ''%-2 - i*«%-3 * - >«7yA-5 * i’i6yA-6 - 

* 1387y-.g - 715y'_, ^ 286y^.^g - TSy^.^^ . I3y^.„ - y^.j, 

‘^uyk = yk - ^%-i * '‘K-2 ■ ^^K-3 * i»®‘yA-4 - * ®»»®yA-6 - ®«%-7 

* 3003y-.g - 2002y'.g * lOOly^.^g - 364y^.jj * Sly'.^j - My^.^, + y^.j^ 

^isJA = yk - ^K-1 * ^^%-2 ■ *^K-3 * i“®>A- 4 - ®““®yA-s * ®“®®y; -s - ««=yA-7 

.6435y^.g - SOOSy^.g . 3003y;.^g - 1365y_;.,j ♦ 455y^.jg - lOSy'.jg ^ ISy^.j^ 

■ yA-is 

^le^A = yk -^%-l * l^Oy’.j - seoy^.j ♦ 1820y-.^ - 4368y^_5 * 8008y^.g - 11440y;_, 

* I2870y-., - 11440y^.g . 8008y’.jg - 4368y^.^^ . 1820y^_jg - seOy'.jj 
-120yi.l4-16y^.l5-y^.l6 

^liyk = yk - * ^^K-2 - ®«“yA-3 * - «l««yA-5 * “376y'.g - 19448y^., 

+ 24310y'.g - 24310y^.g * 19448y;_jg - 12376y;.jj + 6188y;.,j - 2380y^.j3 

. 680y^.j, - 136y^.jg ♦ 17y^.jg - y^.^, 


D-2 


= y; - * “^A-2 - ““yn-3 " ’ «««yA-5 " l“««yi-6 

- 31824y^.j + 43758y_;.j - 48820y^.g * «^*®yA-10 ‘ ^ ‘***^^1^-12 

- " ^o^K-u ■ »^%-is * i«yA-i6 - ‘*yn-i7 " y^-w 

\»>n = yA - ‘®yn-i * ”‘yi-2 ■ »”y,;-3 * ’*’*yA-4 ■ i‘«**yA-5 * 

* M582y^_, - 92378y^.g * »3”8yA-lO ' ”582y^-ll ^ 50”%-i2 ' ”“2y^.„ 

* “«%-i4 ■ " ««»yA-i6 - ”iyA-i7 * i»yA-i8 - yA-w 

^2oyA = yA - =*“yA-i * ^*®yA-2 - “^®yA-3 * ««yA-4 - i«®<yA-5 ^ 3"«»yA-6 

- 7’520yA-7 + 135970yj;.j - 167960y^.j + 184756yj;.j5 - 167960y^_jj 

* 135»»0y’.ij - 77520y^.j3 + 38780y^.j^ - 15504y^.j5 * ««VA-16 

-1140y^.l7.190y’.j8-20y^.l9*y^.2o 


NASA-JSC 


D-3 


